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The complex Langevin method aims at performing path integral with a complex action numerically 
based on complexification of the original real dynamical variables. One of the poorly understood 
issues concerns occasional failure in the presence of logarithmic singularities in the action, which 
appear, for instance, from the fermion determinant in finite density QCD. We point out that the 
failure should be attributed to the breakdown of the relation between the complex weight that 
satisfies the Fokker-Planck equation and the probability distribution associated with the stochastic 
process. In fact, this problem can occur, in general, when the stochastic process involves a singular 
drift term. We show, however, in a simple example that there exists a parameter region in which 
the method works, although the standard reweighting method is hardly applicable. 
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I. INTRODUCTION 

The path integral formulation plays an important role 
in nonperturbative studies of quantum field theories due 
to the possibility of using Monte Carlo methods. The 
basic idea is to generate field configurations with a prob¬ 
ability e~^ and to evaluate the path integral by taking a 
statistical average. It is not straightforward, however, to 
apply such an approach to cases with a complex action 
S since one can no longer view e“'® as the probability 
distribution. This “complex action problem” occurs, for 
instance, in QCD at finite density or with a theta term, 
in Chern-Simons gauge theories, and in chiral gauge the¬ 
ories. It also appears in supersymmetric gauge theories 
and matrix models, which are relevant in nonperturba¬ 
tive studies of superstring theory. 

The complex Langevin method (CLM) Hi attempts 
to solve this moblem by extending the idea of stochastic 
quantization [3| for ordinary systems with a real action to 
the case with a complex action. This necessarily requires 
complexifying the real dynamical variables that appear 
in the original path integral. A stochastic process for the 
complexified variables is defined by the Langevin equa¬ 
tion with the complex action, and expectation values in 
the original path integral are calculated from an average 
of corresMnding quantities over the stochastic process. 
See Ref. Q for a pedagogical review on this method. 

One of the recent developments in this method is the 
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clarification of a necessary condition for convergence to 
correct results Q. In fact, there is a lot of freedom in 
defining the stochastic process that corresponds formally 
to the original path integral. By using this freedom, one 
can try to satisfy the conditions for convergence. The so- 
called gauge cooling @ is a new technique of this kind, 
which has made it possible to study finite density QCD 
in the high temperature (deconfined) phase 0. 

Despite these developments, there is still a puzzle con¬ 
cerning the treatment of logarithmic singularities in the 
action For instance, the effective action for QCD 

involves the logarithm of a determinant, which represents 
the effects of dynamical quarks. At finite density, the de¬ 
terminant A becomes complex in general, which causes 
the complex action problem. In this case, the effective 
action has an ambiguity due to the branch cut of the 
logarithm; the phase of the determinant can be defined 
only modulo 27r. For instance, one may use the drift 
term cllogA = A“^cIA in the complex Langevin equa¬ 
tion, which corresponds to regarding the — log A term in 
the action as a multi-valued function. It has been found, 
however, that the CLM with this prescription can give 
wrong results in a simplified model when the phase of 
the determinant rotates frequently during the stochas¬ 
tic process [1|. The wrong results turn out to be close 
to (but not equal to) the results obtained for the phase- 
quenched model, in which the phase of the determinant is 
quenched. This problem has not shown up yet in recent 
QCD simulations at finite density PH- 

Here we provide new insights into this issue from the 
viewpoint of the Fokker-Planck (FP) equation. We start 
with simple examples and discuss more general cases to¬ 
wards the end. 
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II. ONE-VARIABLE CASE 


evolution of p follows the usual FP equation 


A. Formulation of the CLM 

Let us consider a simple example defined by the parti¬ 
tion function 


d 

-^P = Lop , ( 8 ) 

= -k + i) ' 


Z = J dxw{x) , w{x) = {x + iaY e “ 


( 1 ) 


where a; is a real variable and a and p are real parameters. 
For a ^ 0 and p ^ 0, the weight w{x) is complex, and 
the idea of important sampling cannot be applied to o 
by regarding w{x) as the Boltzmann weight. 

Following the usual procedure in CLM, we define the 
drift term 


v{x) 


w{x) ^ 


dw{x) 

dx 


P 

X + ia 


- X , 


( 2 ) 


and complexify the variable as x z = x + iy. The 
action 


S{z) = — logu;(z) = —p\og{z -\- ia) -\- z^ (3) 


after the complexification has a logarithmic singularity 
at z = —ia for p ^ 0, which causes the aforementioned 
ambiguity due to the branch cut. However, we empha¬ 
size that this is not an issue, in general, because all we 
need in formulating the CLM, as we show below, is the 
single-valuedness of the drift term v{z) after the complex¬ 
ification and the single-valuedness of the complex weight 
w{x) as a function of x. These are satisfied in the present 
case even for a non-integer p. 

The complex Langevin equation corresponding to the 
partition function © can be written as 


^ = v{z) + r]{t) = - z + pit) , (4) 

at z + la 

where p{t) represents a real Gaussian noise satisfying 
{p{t)p{t')) = 2 6it — t'). We dehne the probability dis¬ 
tribution P{x, y; t) of x{t) and y{t) at the Langevin time 
t. Its time evolution follows the FP-like equation 



LP 

— |-(Re^;W,, + —|p 
d 

+ -^{-{ixav)^^:^+iyP} . 


( 5 ) 

( 6 ) 


The crucial point in the CLM is that there exists a 
complex weight p(x; t), which is related to the probability 
distribution Pix, y; t) through 


J 0{x)p{x;t)dx = j 0{x + iy)P{x,y;t)dxdy (7) 

under certain conditions where 0{x) are observables 
that admit holomorphic extension to 0{x + iy). The 


which has a time-independent solution pix\t) oc wix), 
with wix) given in o since it is annihilated by the op¬ 
erator in parenthesis in ([9]). Thus, the necessary and 
sufficient conditions for being able to calculate the ex¬ 
pectation value of O with respect to ([T]) by the CLM 
are: 

i) The relation ([7|) between p and P holds. 

ii) The solution p{x-, t) of the FP equation ® asymp¬ 
totes to wix) as t —>■ cxD up to some constant factor. 

As possible observables in the present example o , we 
consider O = x^, where A: is a positive integer. Assuming 
the ergodicity of the stochastic process, the right-hand 
side of 0 at t = oo can be evaluated by taking the 
time average of z(t)^, where z(t) is obtained by solving 
©■ We find numerically that this method gives correct 
results only for sufficiently large |a| for each p. In what 
follows, we clarify the reason why it fails at small |a|. 


B. Spectrum of the “FP Hamiltonian” 

First we have investigated numerically the eigenvalue 
spectrum of the “FP Hamiltonian” (—Lq) defined by ([9]), 
assuming that the complex weight pix) falls off rapidly 
as |a;| —>■ oo. 

As is clear from what we wrote above, we have an 
eigenfunction pix) = wix) with zero eigenvalue for ar¬ 
bitrary p and a. When p is a positive odd integer and 
a = 0, we have another zero mode pix) = 

For p > I, negative eigenvalues appear in the small |a| 
region. (Note that, when a = 0, we have an eigenfunc¬ 
tion pix) = which corresponds to the smallest 

eigenvalue A = —(p— I) for any p.) Thus, we find that the 
desired solution pix\t) oc wix) is obtained in the long¬ 
time limit of the FP equation (|S]) at arbitrary a for p < I 
and at sufficiently large |q;| for p > I. 

In the parameter region where we have negative modes, 
the complex weight diverges as pix\t) oc 
where A„iin is the smallest eigenvalue and Pmin(a;) is 
the corresponding eigenfunction. Clearly this behavior 
is incompatible with the relation 0 considering that 
Pix,y;t) > 0 and / dxdy Pix,y;t) = 1. This implies 
that the relation 0 between p and P must be broken at 
least in this region. 

Applying this kind of argument to a general multi- 
variable case, we find that the condition like ii) in section 
III Al is automatically satisfied if the condition like i) holds, 
as far as the probability distribution P asymptotes to a 
unique function. 
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C. The relation between p and P 

Let us then consider what can to wrong with ([7]). In 
the derivation of © given in Ref. the authors use 

J 0{x + iy)P{x,y;t)dxdy 

= jo(. + iy,m..v,0)dxiy. ( 10 ) 

where 0{x + iy\t) is defined by solving 


let us take the nth derivative of the above expressions 
with respect to t at t = 0. Using (HU, we obtain terms 
such as 


lim 

T —>-0 . 


{L^YO{z) 

jz + iaj^ 


P(x,y;t)dy , 


(16) 


which diverges as for n > 1 if P{x, y; t) is 

nonzero at (x^y) = (0, —a). 

In order to describe the actual situation, let us define 
the radial distribution 


^0{x + iy-t) = L^O{x + iy;t) (11) 

with the initial condition 0{x + iy; 0) = 0{x + iy). The 
symbol represents an operator satisfying {L^ f,g) = 
{f,Lg), where {f,g) = f f{x,y)g{x,y)dxdy, assuming 
that / and g are regular functions with sufficiently fast 
fall-off as 1x1, |j/| —5> oo. 

In order to prove (dni), they consider [I 

-F('r) = y 0{x + iy;T)P{x,y;t-T)dxdy (12) 


for 0 < T < t, which interpolates both sides of Eq. (unj. 
Taking the derivative with respect to r, they get 

^-F('r) = J ^L^O{x + iy; r)|p(x, y; t - T)dxdy 

—J 0{x + iy;T)LP{x,y;t — T)dxdy . (13) 


Naively, the two terms cancel through integrating by 
parts, which implies that F(t) is independent of r and 
hence (fTOl) . In order to justify the partial integration, 
however, one should be able to neglect the boundary 
terms. This requires that P{x, y; t) decreases sufficiently 
fast as |x|, |y| —>• oo Q. In addition to this requirement, 
one also needs the holomorphy of the drift and of the 
observables to prove ©• 


D. Diverging boundary terms due to singularity 


In the present example, the fast fall-off of P{x,y;t) 
as |x|,|j/| —>■ oo is satisfied due to the — z term in 
0- However, we should be careful of the singularity 
at {x,y) = (0, —a). In order for the boundary terms to 
be neglected, it is required that the limits 


lim 

y^-a 


lim 


0{z;t) 


P{x,y;t 


iy + a)!' yy’.y^ H(x, y; t 
J \z + ia\^ 


T)dy\^ , 
— T)dx 


(14) 

(15) 


should exist for arbitrary t and r. Note, in particular, 
that 0{z; t) obtained by solving (fTT|l is highly singular 
at z = —ia since the operator L appearing on the right- 
hand side of (HU) involves the singularity. For instance, 


'^(^) = TT— oo) S{\/x‘^ + (y + - r)dxdy (17) 

ZTrr J 

around the singular point {x, y) = (0, —a). For small |q)|, 
we observe that (p{r) ^ r at small r as long as p is not very 
large. In this case, HU still diverges for sufficiently large 
n, and the relation © between p and P can be violated. 
Indeed we find that the CLM yields wrong results in such 
a parameter region. 


E. Results for large p 

In the partition function ©, it is the prefactor {x+ia)P 
that causes the complex action problem. In view of this, 
one might think that the CLM fails when the phase of 
(z(f) -I- ia)P rotates frequently during the time evolution 
from ©• We find that this is not necessarily the case. 

In order to demonstrate this point, we present our re¬ 
sults for large p. Figure 1 shows that the CLM reproduces 
the exact results for |a| > 14 at p = 50. From Fig. 2, 
we find for a = 14 that (p{r) = 0 at r < 6, although we 
observe that the phase of (z(f) -I- *a)®° rotates frequently 
during the stochastic process. 


p=50 



FIG. 1: The real part of the expectation value (x^) obtained 
by the CLM (|4]) is plotted against a for p = 50. The solid 
line represents the exact result obtained analytically by the 
Gaussian integration. 
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p=50 



other hand, from the studies of the eigenvalue spectrum 
of {—Lq), we find that p{x]t) oc is obtained in the 

long-time limit of the FP equation (jS]) for arbitrary a 
with /3 = ±1. Therefore, the failure of the CLM at small 
|a| should be attributed to the violation of © due to 
the singularity in the drift term v{z) = 2j3{z + ia)~^ — z. 
This is also confirmed from the behavior of the radial 
distribution (HZl). 

III. TWO-VARIABLE CASE 

Our argument applies also to the case with multiple 
variables. To make this clear, let us consider a case with 
two variables given by 


FIG. 2: The radial distribution defined in (I17II is plotted 
for p = 50 with a = 13 (solid line) and a = 14 (dashed line). 


We would also like to mention that the complex action 
problem is extremely severe at a = 14 and p = 50. As 
a standard quantity that measures the severeness of the 
complex action problem, let us consider 


R = 


w{x) \ 
\w{x)\/ Q 


{x + iaY\ Z 

|a; -I- ia\P )^ Zq 


(18) 


where Zq is the partition function of the phase-quenched 
model Zq = J dx |rc(a;)| and the expectation value ( • )o 
is taken with respect to it. In the present case, both Z 
and Zq can be calculated analytically by performing the 
Gaussian integration. We find that R ~ — 7.4 x 10 “® 
at a = 14 and p = 50 . One can imagine how hard 
it is to obtain correct results if one performs a Monte 
Carlo simulation of the phase-quenched model and ap¬ 
plies the standard reweighting formula to obtain the ex¬ 
pectation values with respect to the original partition 
function ©• Thus, the advantage of the CLM over the 
reweighting method can be appreciated even in this sim¬ 
ple one-variable case. 


F. Non-logarithmic singularities 


Another interesting implication of our argument is that 
the possible failure of the CLM is not specific to logarith¬ 
mic singularities in the action. Indeed we have found that 
the CLM can fail for the weight w{x) = with the 

action 


S{x)=P{x + ia) ^-l-x^/2. (19) 

Note that the action S{z) after complexification does not 
involve a logarithmic singularity, which means, in partic¬ 
ular, that there is no issue of ambiguity associated with 
the branch cut. Nevertheless, we find that the CLM fails 
s-t l®l ^ 1-2 for /? = 1 and at |a| < 1.7 for /3 = —1. On the 


Z = J dxidx2w{xi,X2) , (20) 

W{XI,X2) = (Xi +lX2)^’e-("^)'/2-(.2-a)V2 ^ (21) 

where xi and X2 are real variables. The parameter a is 
real, while p is a positive integer. 

We have studied numerically the eigenvalue spectrum 
of the operator (—Lo)i assuming that the complex weight 
p{x\,X 2 ) falls off rapidly as {xiY + ( 2 : 2 )^ —t 00 . First we 
obtain the desired zero mode p{xi,X 2 ) = w(xi,X 2 ) for ar¬ 
bitrary p and a. When a = 0, we have another zero mode 
p = \xi-\-ix 2 \^ / 2 -( 2 : 2 ) /2 for any p. Forp>l, neg¬ 

ative modes appear at small |a|. (Note that, when a = 0, 
we have an eigenfunction p = {xi -\-ix2) / 2 -(a: 2 ) / 2 ^ 

which corresponds to the smallest eigenvalue A = —(p—1) 
for any p.) Thus, we can make an argument analogous to 
the one-variable case ©. Indeed we find for p = 1,2, 3 
that the CLM with complexified variables zi and Z 2 gives 
wrong results at small |a|. This can be understood from 
the behavior of the radial distribution for r = \zi + iz 2 \- 


IV. IMPLICATIONS TO FINITE DENSITY QCD 

Let us discuss the implication of our argument to finite 
density QCD, which involves the complex fermion deter¬ 
minant det(Zl -|- m) in the partition function, where D 
represents the Dirac operator and m is the quark mass. 
The determinant can be written as the product of the 
eigenvalues Afe of {D -\- m). The drift term of the com¬ 
plex Langevin equation involves where d 

represents the derivative with respect to the complexified 
gauge field. 

According to our argument, the problem we discussed 
does not appear as long as the distribution of Afe is prac¬ 
tically zero at the origin, even if the phase of the fermion 
determinant rotates frequently during the stochastic pro¬ 
cess. (See our results for large p in the one-variable case.) 
This is consistent with the results of recent QCD simu¬ 
lations at finite density, where the distribution of Afe has 
the desired property due either to large quark mass m 
or to high temperature [l3| . 
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On the other hand, the eigenvalues of D obtained in the 
CLM are speculated to accumulate at the origin in the 
chiral limit (corresponding to the m ^ 0 lirnit) when the 
chiral symmetry is spontaneously broken |13l| . If true, 
the CLM will have problems in that parameter regime 
unless some new idea is invoked. 


V. SUMMARY 

We have discussed the issue in the CLM concerning 
the logarithmic singularities in the action. The standard 
drift term corresponds to regarding the logarithm in the 
action as a multi-valued function of the complexified vari¬ 
ables. The CLM with this drift term is known to give 
wrong results in some cases. Theoretical understanding 
of this problem is important, for instance, in applying 
the method to finite density QCD in the low tempera¬ 
ture (confined) phase with light quark mass. 

First we emphasized that the multi-valuedness of log¬ 
arithmic terms in the action cannot be considered the 
cause of the problem since one can formulate the method 
without referring to the action, as we have done in sec¬ 
tioning This is also indicated by the example in section 

EB 

Rather, the problem should be attributed to the possi¬ 


ble breakdown of the key relation between p and P due to 
the singularities in the drift term of the complex Langevin 
equation. In particular, we pointed out that the relation 
can be violated due to the boundary terms appearing 
from integrating by parts in proving OT, which diverges 
unless P is practically zero around the singularities. This 
assertion was supported by simple examples. A more 
quantitative analysis will be reported in the forthcoming 
publication. 

The “FP Hamiltonian” can have negative modes only 
if the key relation between p and P is violated. Note, 
however, that the key relation can be violated even if the 
“FP Hamiltonian” does not have negative modes. Hence, 
the appearance of the negative modes should be regarded 
merely as an indicator of the violation of the key relation, 
the latter being the cause of the problem. 

To conclude, we hope that the new insights gained in 
this work will be useful in developing the method further 
in cases with singularities in the drift term. 
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